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Abstract 

After presenting a new approach to separate the total angular momentum of an electromagnetic 
beam into the spin and orbital parts, which manifests that the spin angular momentum originates 
from that part of the linear momentum density the total amount of which is equal to zero, I 
show that the orbital angular momentum as well as the spin angular momentum of a non-paraxial 
monochromatic beam is dependent on the polarization ellipticity a. The cr-dependent term of the 
orbital angular momentum is mediated by the recently advanced symmetry axis I; and the trans- 
verse component of the orbital angular momentum is consistent with the transverse displacement 
of the beam's barycenter from the plane formed by I and the propagation axis. For a beam of 
angular-spectrum scalar amplitude f(k p ,ip) = fo(k p ) exp(il(p), the total angular momentum per 
unit energy in the propagation direction is simply for = ^ and ~ for = 0, where I is an 
integer, to is the angular frequency, and is the angle between I and the propagation axis. 

PACS numbers: 42.60.Jf, 42.25.-p 
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I. INTRODUCTION 



It is well established that electromagnetic beams can carry angular momenta. It was 
shown [l[ [2] that the total angular momentum (TAM) could be separated into the spin and 
orbital parts. The spin angular momentum (SAM) 31 and the orbital angular momentum 

n fl ft n n 

B)AM) [J, have been experimentally measured. The transfer of SAM |6j, [3, ^| and OAM 
been observed. The distinction between SAM and OAM in the 



14 



id, 



13|. 



interaction with particles has also been demonstrated [111 . [12j. The classical mechanics 
aspect of the angular momentum opens the door to the optical micromachines 
The quantum mechanics aspect opens the door to the quantum information 

It is commonly believed on the basis of the knowledge of paraxial Laguerre-Gaussian 
beams [l| that the SAM and the OAM of free beams are independent of each other and 
are carried 2], 5), respectively, by the polarization and the helical phase factor exp(i/0), 
where I is an integer. But a recent experiment [l7] converted partly the SAM into the OAM 
by a high numerical aperture. This fact might signify some connection between the SAM 



and the OAM 



181 ] . The purpose of this Letter is to investigate this connection for non- 



paraxial vector beams. The non-paraxial beam that I examine is characterized by the newly 

n 

advanced symmetry axis I [19] that makes an angle O with the propagation axis. The case 



of© 



20 



21| represents the uniformly polarized beam (including the Laguerre-Gaussian 



beam) in the paraxial approximation. The case of O = [22] represents the cylindrical 



vector beam 
transformed 



23 



25 



24j. And the beam of O that is neither equal to | nor equal to was 
261 ] from a beam of O = | by the transmission at an interface between two 
different dielectric media. It will be shown that both the SAM and the OAM are dependent 
on the polarization. The polarization-dependent term of the OAM is mediated by I. It will 



also be shown that t 
displacement 



19 



re transverse component of the OAM is consistent with the transverse 



25 



261 ] of the beam's barycenter away from the plane formed by I and 



the propagation axis. To this end, I will first show in a new attempt to separate the TAM 
that the SAM originates from that part of the linear momentum density the total amount 
of which is equal to zero. It is this property that makes the SAM independent of the choice 
of the origin. 
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II. SEPARATION OF THE TAM 



Consider an arbitrary electromagnetic beam in free space. Its electric vector in the 
position space can be written as the following integral over the plane wave, 



where 



( k \ 



(1) 



yk zJ 

is the wave vector each element of which is real, Eq^ouj 2 = k 2 , and E(k) is the electric 
vector in the wave-vector space. Obviously, integral expression ([T]) leads to the following 
t r ansf or mat ion , 

c(-k) = -w(k), E(-k) = E*(k), (2) 

where the superscript * denotes the complex conjugate. The magnetic vector of the beam 
is derived from Eq. ([T]) and the Maxwell equation to be 



H(x,t) = - 



1 



k x E 



exp[i(k • x — Lut)]d 3 k + c.c. > . 



(3) 



2 [ (2^)3/2^0 J oo 

According to Eqs. ([1]) and Q, the linear momentum density p = e fi £ x 7i is separated 
into two parts, p = Pi + P2, where 

£0 f E' • E 



Pi 



P2 



4(2tt) 3 
4(2tt) ; 
4(2tt) j 
4(2vr) ; 



UJ 

E' • E* 



_ ke i(k'-k)-x e -i(u,'-u,)t d 3 k , d 3 k + c _ c 



UJ 



(4) 



E'-k 

uj 
E'-k 



Yie i{M'+k)^ e -i{u 1 '+u 1 )t d 3 kl(i 3 k 



E * e K^-^ e -i^'-^t d 3 k f d 3 k + 



UJ 



(5) 



E = E(k), E' = E(k'), oj = uj(k), and uj' = uj(k'). Correspondingly, the TAM 3 = J xxpd 3 x 
is also separated into two parts, J = L + S, where j^] 



L 

S 



x x p^ 3 x 
x x p 2 d 3 x 



£0 



E f (k x V k )Ed 3 k, 



— E* x Ed 3 k, 

iuj 



(6) 
(7) 



+ e. 



+ e 2 



and the superscript f stands for the conjugate transpose. The 



details to derive the above equations are given in Appendix. Eq. ([7]) shows that S arises 
from that part of the linear momentum density the total amount of which is equal to zero, 



/ 



p 2 d 3 x = 0. 



(8) 



Consequently, the spatial translation x — ► x — x does not change S. It is thus tempting 
to regard S and L as the SAM and the OAM, respectively. Property (jHD also tells us that 



there is no paradox [28|, |29j about the SAM of the plane wave. 



III. ANGULAR MOMENTUM OF MONOCHROMATIC BEAMS 

The electric vector of a monochromatic wave propagating in positive z direction can be 
expressed in terms of the angular spectrum as 

£(x,t) = ^j^- J E(k x , k y ) exp[i(k • x - ut)]dk x dk y + c.c. j , (9) 

where k 2 = k 2 — k 2 — ky, and the integration limit k 2 + k 2 < k 2 is omitted for brevity. 
Similarly, one has the following transformation from the integral expression (Q, 

kz ( k x , ky) k z (yk x , kyj, 

u{-k x ,-k y ) = -u(k x ,k y ), (10) 
E( k x , ky) = E {k x , ky), 

which is consistent with transformation In much the same way as before, one finds that 
the linear density of the TAM, the integral of the angular momentum density over the cross 
section z = constant, is independent of the z coordinate and is separated into the orbital 
and spin parts in a natural way, J = j x x pdxdy = L + S, where 



£o 



, _ dE , _ dE / ^ dE , ^ dE 
-e x k z E* ■ — + e y k z E* ■ — + e A k x E* ■ — - k y E ■ — 
oky ok x \ ok y ok x 



S = — I E* x Edk x dk y . (11) 



It is easy to show that the above equations can be obtained directly from Eqs. (jSj) and <^ 

dE 



by noticing -jP- = and replacing the triple integral with the double integral. 
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IV. PROPERTIES OF THE ANGULAR MOMENTUM OF NON-PARAXIAL 
BEAMS 



A representation theory for non-paraxial vector beams was recently advanced in Ref. [19 ]. 
The electric vector in the position space is given by Eq. The electric vector E of the 
angular spectrum is factorized into three factors, 

E = maf, (12) 

where m is the 3x2 mapping matrix specified by a newly found symmetry axis I, 

'-(:) 

is the Jones vector satisfying the normalization condition a^a = 1 and describing the polar- 
ization state of the angular spectrum through the polarization ellipticity, 

a = oraa, = — i(a*oi2 — a^cti), 

(, ,\ 

a = is the Pauli matrix, and / is the scalar amplitude of the angular spectrum. 

V / 

Letting I lie in the plane zox and make an angle G with the propagation axis, 

I = e z cos G + sin 0, 



the mapping matrix takes the following form [19], 



1 

m 



' (ky + kf) sin 6 — k z k x cos 6 kk y cos© ^ 



fc|k x I] 



(14) 



—k y (k z cos + k x sin Q) k{k z sin O — k x cos G) 
, (k% + ky) cos O — k z k x sin G — kk y sinQ j 

where |k x I] = [k 2 — (k z cos G + k x sin G) 2 ] 1 / 2 . The a is assumed to be the same to all the 
elements of the angular spectrum. In order that the Laguerre-Gaussian beams be included, 
the scalar amplitude / is chosen to have the following form, 

f(k P , <fi) = fo(k p ) exp(^), (15) 

in the circular cylindrical system in which k = k p + k z e z , k p = k p e p = k x e x + k y e y , 
k x = k p cosip, k y = kpSiiiLp, k z = [k 2 — k 2 ) 1 / 2 , e p and e 9 are, respectively, the unit vectors 
in the radial and azimuthal directions, where fo(k p ) is square integrable. 



Let us first look at the property of the SAM component in the propagation direction, 

S z = ^ f (E*E y - E* y E x )k p dk p d V . 
Substituting Eqs. (|T2|) - (|T5|) . one has 

S z = 2ne - I hf (k p )\%dk p , (16) 
to J k 

which does not depend on the beam parameter 0. This expression has the validity that does 
not rely on the paraxial approximation. In the linear paraxial approximation in which fo(k p ) 
is sharply peaked at k p = and therefore the factor k z in the integrand is approximated as 
k z « k, one arrives at 

S z « -W, (17) 

n 

which was first observed in Ref. [lj, where 

W = 2ne [ \f (k p )\ 2 k p dk p (18) 
Jo 

is the linear density of the energy, the integral of the energy density over the cross section. 
Similar considerations lead to 

S x = S y = (19) 

for the SAM components in the x and y directions. 

Then let us investigate the OAM component in the propagation direction, 

L z = — f ^{-i-^-)m p dk p d V . 

Substituting Eqs. (fT2l)- ([15l) and performing a lengthy calculation, one finds 



i -> +2 ^f (K l+ ^S) 4}iwv*„ ( 2 o) 

where |0| < | is assumed 0]. This expression also has the validity that does not rely on 
the paraxial approximation. It shows an apparent dependence on the beam parameter B 
as well as the polarization ellipticity a. If G = | which describes the uniformly polarized 
beams (including the Laguerre- Gaussian beams) in the paraxial approximation, one has 

L z = l ~W + 2ne ^ j* (l - j^j \fo(k p )\ 2 k p dk p . (21) 
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In the linear paraxial approximation, it reduces to 

I 



UJ 



W. (22) 



This re 

beams 



ation was also obtained in Ref. If @ = which describes the cylindrical vector 



2j, Eq. (J2DD becomes 

L z = -W- 2ne - [ hf {k p )\ 2 k p dk p . (23) 

UJ UJ J K 

In the linear paraxial approximation, it reduces to 

L, « '-^W, 

UJ 

which is cr dependent. 

Eqs. ( Tl6|) and ( 120|) indicate that the TAM component in the propagation direction is 

J, = -W + toe Q - f \ ( 1 + ^'^ni ) \fo(K)\%dk p . (24) 
uj uj J 2 \ \k z — kcosB|/ 



If 6 = |, it reduces to 



If = 0, one simply has 



J z = —W. (25) 

UJ 



J z = -W, (26) 

UJ 



which is cr independent. 

With the above conclusions, it is ready to explain the spin-to-orbital angular momentum 



conversion through focusing a beam of 9 = | [171 ] . Eq. fl23|) shows that the TAM per unit 
energy of this beam in the propagation direction is Before focusing, the SAM and OAM 
per unit energy of the paraxial beam in the propagation direction are approximately - and 
-j, respectively, as Eqs. ffTTl) and fl22|) show. After focusing, the SAM per unit energy of the 
non-paraxial beam is obtained from Eq. (TlBI) to be 

cr jQ(k z /k)\fo(k p )\ 2 k p dk p 



u f k \fo(k p )\%dk p 

indicating that only a fraction of the incident SAM remains in the focused beam, where 
fo(k P ) of the focused beam is not sharply peaked so that the linear paraxial approximation 
is no longer applicable. Since a lossless focusing system cannot change the TAM, the rest of 
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the incident SAM is converted into the OAM of the focused beam [30J. In fact, the OAM 
per unit energy of the focused beam is obtained from Eq. (|2ip to be 

rk, 



I a 
- + - 

uo uo 



1 - 



f K (k z /k)\f (k p )\%dk p 

Jo \fo{k P )\ 2 k p dk p 



At last, let us discuss the OAM component in the x direction, 



£o 



d 

k z W(i—)Ek p dk p dip. 



(27) 



Substituting Eqs. (fl2]) - (fl5]) and performing a similar calculation, one gets 



-2ne - 



a cot O 



~ 1 + 



k z — k cos O 
\k, — k cos 01 



\fo(kp)\ 2 kpdk p , 



(28) 



which is also dependent on as well as a. This is surprisingly contrary to the common 
belief 17] that the OAM is carried only by the helical phase factor. If = 0, one has 



L x = 0. On the other 
one approximately has 



land, if |0| ^> 56 with 56 being the divergence angle of the beam, 
o cot , 



19J 



-W. 



UJ 



As a matter of fact, it can be seen from Eqs. ( 1271) and ( fl8i) that in the linear paraxial 
approximation, the OAM per unit energy in the x direction, is simply related to the 
transverse displacement [3] of the beam's barycenter from the plane formed by the 
symmetry axis I and the propagation axis in the following way, 

w - & < 29 > 

A similar calculation shows that the OAM component in the y direction vanishes, 

L y = 0, (30) 

indicating together with Eqs. (1201) and ([25]) that the OAM is located in the plane formed 
by the symmetry axis I and the propagation axis. 

V. CONCLUDING REMARKS 



In summary, I presented a new approach to the separation of the TAM into the spin and 
orbital parts and found that the spin originates from such a part of the linear momentum 
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density the total amount of which vanishes. Contrary to the common belief, I showed that 
both the OAM and the SAM are dependent on the polarization ellipticity a. Their difference 
in respect of the cr-dependence lies in the different effects of the symmetry axis I. The a- 
dependent term of the OAM is mediated by I as Eqs. (12"U1) . f[2"8"j) . and f[3~0~j) indicate, whereas 
the SAM is free of the effect of I as Eqs. (I16p and (TiT?|) indicate. The mediation of the 
a-dependent term of the OAM by I may offer further insights into the nature of the angular 
momentum of light. 
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APPENDIX: DETAILS OF THE DERIVATION OF EQS. © AND © 

Let us first derive Eq. fl§]). Substituting Eq. (j3J into L = Jx x pid 3 x, one has 



L = Li + L 2 + c.c, 



(A.l) 



where 




(A.2) 



and 




(A.3) 
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Upon integrating Eq. (1A.2I) over the configuration space and noticing the following properties 
of Dirac's 5 function and its first-order derivative, 

POO POO 

S(t) — — I exp(iut)du, —i5'(t) = — I uo exp(iujt)duj, (A. 4) 

27r J-oc 27r J-oo 

one obtains 

L : = g y (J^e* - fc 2 e y )^-^e^ ( "' +w) *5'(A;; + fe a )5(ifej + k y )8{k' z + k z )d 3 k'd 3 k 
+ f. J (k z e x - k x e z )^-^e~ l ^ + ^8(k' x + k x )8'(k' y + k y )8(k' z + k z )d 3 k'd 3 k 
+ % j ( kxe v ~ K^)^-^e- i{uJ ' +uj)t 5{k' x + k x )5{k' y + k y )5'{k' z + k z )d 3 k'd 3 k. 
It is changed by eliminating the 8 functions into 

U = j. j (k y e z - k z e y ) E{K > ~ h) ' V ^'-^-^ + ^'(£4 + k x )dk' x d 3 k 

+ g J (k z e x - k x e z ) n ' K ^ *v ~ K) ' ^ e -^-^-^5\k' y + k y )dk' y d 3 k 
Noticing the following property of the derivative of the 5 function, 

eta 



/ /(*)$'(* - t Q )dt = -f(t ), *i < t < t 2 , 



(A.5) 



and taking transformation (j2J) into account, the above equation is reduced to 
Li = £o f k y e z - k z e y / gE* + .J^. \ ^ 

Ai J uj V £o/-*o^ / 

+ £o W E -^ + ^E*-eVa; 

,:>) f V(kxV t )E*A, 

u; 

where the superscript T denotes the transpose. By making the variable replacement k — ► — k, 
it is changed into a familiar form, 

L x = - [ — E f (k x V k )Bd 3 k. (A.6) 
4 J iuj 

Since — iV^ is a Hermitian operator, the Li given by Eq. (1A.6P is real. A similar calculation 
produces from Eq. (1A.3j) 

L 2 = Lx. (A.7) 
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It is clear that substituting Eqs. (IA.6[) and (IA.7D into Eq. (lA.ip will yield Eq. f| 
Then we derive Eq. ([7]). Substituting Eq. ([5]) into S = Jxx p2d 3 x, one has 



where 



and 



Si 



4(2vr) ; 



~4(2tt) ; 



S = Si + S 2 + c.c, 



d 3 A;'(i 3 A; / d 3 x— -x x Ee i(k ' +k) - x e" i(w ' +w)t , 



d 3 k'd 3 k I d 3 x— -x x EV^'-^e^'-^*. 
a; 



(A.8) 



(A.9) 



(A.10) 



Upon integrating Eq. (1A.9|) over the configuration space and noticing Eq. (1A.4j) . one obtains 



Si 



Z£0 

4 

ie 



+ — I (E z e x - E x e z 
+ f / (E x e y - E y e x ] 



E'- 


k 


a; 


— 


E'- 


k 


UJ 




E'- 


k 





e -^'+")^'(^ + + k y )S(k' z + k z )d 3 k'd 3 k 

e-^'+^Sik'x + k x )5\k' y + k y )5{k' z + k z )d 3 k'd 3 k 

e -i^'+^)t s ^ + ^(jfc' + + k z )d 3 k'd 3 k. 



It is changed into, by eliminating the 5 functions and taking Eqs. flA.51) and (T5J) into account, 



s ' = f 



z£ /" E y e z - E z e y ^ _ <^ d z k [ Ez&x ~ E ^ z \^ . rf 3 fe 



UJ 



dk r 



4 



UJ 



+ 



4 J w " <% 2 

From the transversality condition k • E* = 0, we know that 



r)E* 

-^■ k '«: = - £ *- 



Si thus takes the following form, 

Si = - f ^E* x Ed 3 k, (A.ll) 
which is clearly real. Similarly, S 2 in Eq. (1A.10I) is found to be real and is equal to Si, 

S 2 = Si. (A.12) 



Substituting Eqs. flATTTj) and flA~T2|) into Eq. (|AT8|) will yield Eq. © 
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